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ABSTRACT 

Gamma Ray Bursts (GRB) observed up to redshifts z > 8 are fascinating ob- 
jects to study due to their still unexplained relativistic outburst mechanisms and 
a possible use to test cosmological models. Our analysis of 77 GRB afterglows 
with known redshifts revealed a physical subsample of long GRBs with canonical 
plateau breaking to power-law light curves with a significant luminosity L* x - break 
time T* correlation in the GRB rest frame. This subsample forms approximately 
the upper envelope of the studied distribution. We have also found a similar rela- 
tion for a small sample of GRB afterglows that belong to the intermediate class 
(IC) between the short and the long ones. It proves that within the full sample 
of afterglows there exist physical subclasses revealed here by tight correlations of 
their afterglow properties. The afterglows with regular ('canonical') light curves 
obey not only a mentioned tight physical scaling, but - for a given T* - the more 
regular progenitor explosions lead to preferentially brighter afterglows. 

Subject headings: gamma-rays bursts: general - radiation mechanisms: nonther- 
mal 
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Introduction 



The detection of GRBs up to high redshifts (z=8.2; ISalvaterra et a 



( 120091 )). larger than Supernovae la (SNela) (z ma: , 
objects appealing for possible use in cosmology. 



1.77; 



Riess et al 



■l (l2009h . lTanvir et al. 
(EooJ), makes these 



The problem is that GRBs seem not to 



20051 : iLiang et all 12 009; 



thermore Cardone et al. 



way, several GRB luminosity indicators ( 


Amati et al. 


2008 


; Fenimore 


& Ramirez - Ruiz 


2000; 


Norris et al. 2000: 


Ghirlanda et al. 


2004 


Liang & Zhang 


2005. 


2006; 


Ghirlanda et al. 


2006) 


and their use to constrain cosmological parameters 


( Firmani 


et al. 


2006a: 


Liang & Zhang 



Qi et al. 1120091 : llzzo et al.l 120091 ) have been proposed till now. Fur 



(j2009l ) have derived an updated GRBs Hubble diagram using the 



l og L y -log T* ('LT') correlation with five other 2D GRBs correlations used by Schaefer 
(200 7J). However, the problem of l arge data scatters in the considered luminosity relations 
( IButler et al.l l2Q09t lYu et al.l 120091) and a possible impa ct of detector thresholds on cosmo- 
l ogical standard cand les (jshahmoradi fc Nemiroff 2009 ) have been di scussed controversially 
( ICabrera et al.ll2009l ). Among these attempts, iDainotti et al.l ( 120081 ) have proposed a way 
to standardiz e GRBs as dist a nce i ndica tor with the disc overy of the LT anti-correlation, 
confirmed by iGhisellini et al.l ( 120091 ) and lYamazakil ( 120091 ) . The fitted po wer-law relation is 
log L* x = log a + b ■ log T*; the constants a and b are determined using the iD'Agostinil (120051 ) 
method. 

In this Letter we study the LT correlation using the extended GRB data set and we 
demonstrate the existence of a physical LT scaling for 'canonical' light curves in the GRB 
rest frame. The regular lightcurve afterglows conform rather tightly to this scaling, while 
the more irregular ones are systematically fainter. A similar correlation is revealed for a 
subsample of GRB afterglows that belong to the intermediate class (IC). Revealing these 
physical correlations can help the (still unclear) interpretation of the physical mechanisms 
responsible for the GRB X-Ray afterglow emission and can infer important information 
about the nature of the emitting source. 



2. Data selection and analysis 

We have analyzed a sample of all afterglows with known redshifts detected by Swift from 
January 2005 up to April 2009, for which the lightcurves include early XRT data and there- 



1 We use the index '*' to indicate quantities measured in the GRB rest frame in which L* x = L* X {T*) is 
an isotropic X-ray lumin osity in the time T* , the transition time separating the afterglow plateau and the 
power-law decay phases (jDainotti et al.l l2008). 



-3 - 



fore can be fitted by a Willingale's phenomenological model ( IWillingale et al.ll2007l ). The 



redshifts z are taken from the Greiner's web page http:/ / www.mpe.mpg.de/'; 



cg/grb.html 



We have compared these redshifts with the values reported by iButler et al.l (120071 ) and we 
find that they agree well apart from two cases of GRB 050801 and 060814, but Butler (pri- 
vate communication) suggested that we should use the Greiner redshifts for thos e two cases. 
For original references providing the redshift data see (IButler et al.l 120071 . 120091 ). Our data 
analysis, includin g derivation of T* and L* x (in units of |s| and [grg/ s] , respectively) for each 
afterglow, follows iDainotti et al.l ( 120081 ) and IWillingale et al.l (120071 ) . The source rest frame 
luminosity in the Swift XRT bandpass, (E min , E max ) = (0.3, 10) keV, is computed from the 
equation: 



, il!n , E max , t) = AirD L (z) F x {E, 



mini E ma r 



where Dl(z) is the GRB luminosity distance for the redshift z, computed assuming a flat 
ACDM cosmological model with Qm = 0.291 and h = 0.697, Fx is the measured X-ray 
energy flux (i n [erg /cm 2 /si) and K is the K-correction for cosmic expansion. Using the 
Bloom et al.l ( 120011 ) expression for the /^-correction and with f(t) being the Swift XRT 
lightcurve we have the relation : 



KFx(E min , E max ,t) = f(t) x 



r 1 
■it 



»/(!+*) 



E${E)dE 



E$(E)dE 



(2) 



where <&(E) is a usual differential photon spectrum assumed to be oc E~ la = E~^ a+1 \ where 
7 a and (3 a are the photon index and the spectral index, respectively. IWillingale et al.l (120071 ) 
proposed a functional form for f(t): 



fit) = f P (t) + fait) (3) 

where the first term accounts for the prompt (the index "p" ) 7 - ray emission and the initial 
X-ray decay, while the second one describes the afterglow (the index "a"). Both components 
are modeled with the same functional form : 



/(*) = 



/•; rxp ( a c - -^j exp (~ ) for t < T< 



for t > T r 



(4) 



FA T, 



exp \ -- 



-4- 



where c = p or a. The transition from the exponential to the power law occurs at the point 
(T c , F c ) where the two functional sections have the same value and gradient. The parameter 
a c is the temporal pow er law decay index an d the time t c is the the initial rise time scale 
(for further details see ( jWillingale et al.l 120071 )). 

For the afterglow part of the light curve we have computed values L* x (eq. [[]) at the 
time T a , which marks the end of the plateau phase and the beginning of the last power law 
decay phase. We have considered the following approximation which takes into accounts the 
functional form, f a , of the afterglow component only: 



/(T a )^/ a (T a )=F a exp 



T 
T 



for t — T„ 



(5) 



where we put the time of initial rise, t a = T p because in most cases the after glow component 
i s fixe d at the transition time of the prompt emission, T p (for details see IWillingale et al. 
(120071 )). Then, with applying eq-s |5] and [2] in eq. [I] one obtains: 



L 



x 



AnD 2 L (z)F x 



(6) 



where Fx = F a exp(— ^f) is the observed flux at the time T a . We have derived a spectral 
index / L for each GRB a fterglow using the Evan's web page |http:/ /www.swift.ac.uk/xrt] 
curves (lEvans et al.l 120091 ) setting a filter time as T a ± <7 Ta ; the T n values togeth e r wit h 
their errorbars, ax a , are derived in the fitting procedure used by IWillingale et al.l ( 120071 ). 
As mentioned above the power law spectrum, <f>(E) oc E~^ a+l \ was fitted with the model 
'phabs*phabs*pow' providing the X-ray spectral index, /3 a . The first absorption component 
is frozen at the Galactic column density value obtained with the NH FTOOL0 and the second 
is the 'zphas' component with the redshift frozen at the value r e porte d in the literature. For 
further details of the spectral fitting procedure see lEvans et al.l ( 120091 ). The lightcurves used 
for the analysis are the same used in lEvans et al.l ( 2009 ). but binned by us in a different way. 



For some of the derived points (L* x , T*) the error bars are large, indicating that the 
canonical lightcurve doesn't fit well the observed lightcurve. We have decided to include 
such cases in the analysis to treat the whole sample in a homogeneous way. Even if points 
with largest - a few orders of magnitude - error bars have no physical meaning, they carry 
information about the light curve irregularity (deviation from the considered model) or 
insufficient amount of observational data for precise fitting. 



2 



http://heasarc.gsfc.nasa.gov/lheasoft/ftools 
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A choice of the Willingale model as a representation for the X-ray GRBs lightcurves 
allows us to use a homogeneous sample of events to study physical correlation in a statistical 
way. Let us point out that the fitting procedure can yield values of Ta within the gap between 
the end of BAT and the beginning of XRT observations, like in the case of GRB 050318. 
One can note that several authors fit the afterglow part of the lightcurve without modelling 
the prompt emission light curve. Thus, they can obtain a nearly perfect power law fits in 
cases where the Willingale model fitting finds a short afterglow plateau phase. 

To analyse how accuracy of fitting the canonical lightcurve, (eq. |3] and @J to the data 
influences the studied correlations we use the respective logarithmic errors bars, at* and 

cry* , to formally define a fit error parameter u = ■^ja\* + <j\* , as measured in the burst rest 
frame. This definition is used to distinguish the canonical shaped light curves from the more 
irregular ones, perturbed by 'secondary' flares and various non-unifo rmities. Th e symm etric 
error bars quoted in the paper are computed with the method of iD'Agostinil ( 120051 ) that 
takes into account the hidden errors and thus gives greater error est imates than the ones 
obtained with the Marquardt Levemberg algorithm (jMarquardtlll963l ). 



Our analyzed sample of 77 GRBs from the redshift range 0.08 — 8.26 includes after- 
glows of 66 long GRBs and 11 GRBs whose nature is debated, the intermediate class (IC) 
between long and short GRBs, described by iNorris et al.l ( 120061 ) as an apparent (sub) class 
of bursts with a short initial pulse followed by an extended low-intensity emission phase. 
Our lon g GRB sample incl udes also 8 X-Ray Flashes (XRFs) ( 060108, 051016B , 050315, 
050319 JGendre et aDl2007h . 050401, 050416A, 060512, 080330 JSakamoto I l200~8h . XRFs 
are scattered within the long GRBs distribution in Fig. [H provi ding further support to a 
hypothesis that both these phenomena have the same progenitors (lloka fc Nakamurall200ll ). 
To study physically homogeneous samples we decided here to analyze the sub-samples of 66 
long GRBs (including XRFs) and of 11 IC ones separately. 



The results 



The obtained l L* x 



versus T*' distributions for long GRBs (Fig. [TJ §) and for a smaller 
sample of IC GRBs (Fig. [2]) clearly demonstrate existence of significant LT correlations, 
characterized in this paper by the Spearman correlation coefficient, p, a non-parametric 
measure of statistical dependence between two variables ( ISpearman 1 1 19041 ) . From a visual 
inspection of Fig. [1] and the analysis discussed later on Fig. [3] one can note that the lowest 



3 see the data table for all long and IC GRBs at 



http://www.oa.uj.edu.p1/M.Dainotti/GRB2010 
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error events concentrate in the upper part of the distribution, forming a highly correlated 
subsample of the full distribution. To visualize this effect we decided to arbitrarily select 8 
points with smallest errors to define our limiting the upper envelope subsample, u < 0.095, 
see the inset panel in Fig. [TJ 

For the full sample of 66 long GRBs one obtains plt = p(logL x ,\ogT*) = —0.68 
and a probability of occurring of such correla t ion by chance within the uncorrelated sample 



P = 7.60 x 10 9 (cf. iBevington fc Robinson I (120031 )). If we remove a few large error points 



by imposing a constraint u < 4 we have a limited sample of 62 long GRBs presented in 
Fig{TJ with plt = —0.76, P = 1.85 x 10~ n and the fitted correlation line parameters^ 
log a = 51.06 ± 1.02 and b = — 1. 06^0 28 > wn il e for the upper envelope sample we obtain, 
respectively, p LT = -0.93, P = 1.7 x 10~ 2 , logo = 51.39 ± 0.90 and b = -1.05^. 

In Fig. [2] we present data for all IC GRBs. One may note that their afterglows are 
characterized with the values of T* in the upper times range of the long GRBs. The IC 



Norris et al. 


2006 


), 060912A ( 


Levan et al. 


2007) 



follow a similar LT relation as the long 
ones. A formally computed correlation coefficient for these 8 GRBs is plt = —0.66. This 
result suggests the existence of an another steeper LT correlation for IC GRBs as the one 
obtained for long GRBs, with different fitted parameters log a = 52.57 ± 1.04 and b = 
— 1.72lo2i- The plotted fit line is below the analogous u < 4 sample of long GRBs showing 
that the IC GRBs and the normal long ones behave differently, but a limited number of 
available IC GRBs inhibits us to make a strong statement in that matter. 

To study the fit error systematic of GRB afterglows we show below, in Fig[3j how the 
limiting upper value for u in the analyzed sample, i.e. how selecting the afterglows with 
increasing precision of L* x and T* fits, influences the LT correlation. We present changes of 
the Plt converging - with decreasing u - toward a nearly linear LT relation, as observed for 
our upper envelope sample. On the figure, e.g., we have 62 long GRBs for u = 4, 33 GRBs 
for u = 0.3, 19 GRBs for u = 0.15, 13 GRBs for u = 0.12 and 8 GRBs left for our limiting 
u = 0.095 . A presented accompanying systematic shift upwards of the fitted correlation - 
as measured in the middle of the distribution as log a — 3.0 ■ b (the fitted correlation line 
at an arbitrarily selected logT^ = 3.0) - with decreasing u, proves that the limiting u C 1 
subsample forms the upper part, the brightest afterglows in the LT distribution. This regular 
trend allows us to conclude that the subclass of all long GRBs with 'canonical' afterglows 



'One may note that the presented fitted slope is different from the slope range quoted in lDainotti et al 



(|2008l ). because in the previous paper in an attempt to reduce the intrinsic scatter in the correlation, the 
authors limited the sample to the GRBs with logL^ > 10 45 erg/s and time parameter log[T* < 5, with a 
possible resulting bias in the fitted correlation parameters. 
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1 1.5 2 2.5 3 3.5 4 4.5 5 



log Ta* 

Fig. 1. — L* x versus T* distribution for the sample of 62 long afterglows with u < 4, with the 
fitted correlation line in black. The upper red line, fitted to the 8 lowest error (red) points, 
forms approximately an upper envelope of the full distribution. The upper envelope points 
with the fitted line are separately presented in an inset panel. 
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2.5 3 3.5 4 4.5 5 5.5 6 

log Ta* 

Fig. 2. — L* x versus T* distribution for the sample of 11 intermediate class GRBs. For the 
picture clarity the three points with very high errors bars (u > 4) are shown without error 
bars. The fit dashed line is presented for the 8 points with indicated error bars, for u < 4. 
Additionally, both fit lines for long GRBs from the Fig. [1] are provided for a reference. 
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•0.75 



a. -0.85 




-0.95 



48.3 



48.2 



48.1 



47.9 



47.8 



47.7 



Fig. 3. — The left vertical axis: correlation coefficients p(logL* x ,\ogT*) versus the error 
parameter u for the long GRBs are presented with black squares. The right vertical axis: 
normalizations of the fitted correlation lines at logT^ = 3.0 versus u are presented with red 
asterisks. 





Fig. 4- 



T* versus the spectral index j3 a for the upper envelope sample, u < 0.095. 
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forms a well defined physical class of sources exhibiting high correl ation of their afterglow 



parameters. Presence of GRBs with light curves deviating from the I Willingale et al.l (120071 ) 
model increases the scatter in the L-T distribution, with larger error points distributed 
preferentially below the small error ones. 

We expect that having a tight LT correlation for canonical light curves, suggesting a 
precise physical scaling to exist between GRBs with different luminosities, should be ac- 
companied with a regular change of the fitted spectral indices (3 a , if any such changes occur 
besides a scatter due to measurement errors. In fact we find the existence of (3 a -\og T* (and, 
of course, (3 a -logL* x ) correlation, see Fig. HI for a sample of small u < 0.095 points, where 
pp aT = p(/3 a logT^) = 0.74, but a non-negligible probability to occur by chance from an 
uncorrelated distribution, P = 1.0 x 10 _1 . It suggests that the (3 a tends to increase for larger 
T*: the fitted relation reads (3 a = 0.43 + 0.19 x \ogT*. For large limiting u = 4 this corre- 
lation becomes weak with ppx ~ 0.16. A warning should be considered. Since we determine 
the value of (3 a with a filter time T a ± ax a the measurement systematics could be introduced 
to the data if (3 varies with time in real objects. To evaluate such possibility we computed 
correlations p((3 a , logT a ) using the observer measured times T a , without a ^-correction, the 
ones used in real f3 a measurements. However - having a large scatter of a values at all T a - 
we do not observe any increase or systematic change of correlations as compared to the ones 
using T* times. This test neither excludes nor supports systematics in the data, leaving the 
possibility that the observed effect could be real. 

To preliminary verify if any redshift systematics exists in the GRB afterglow distribution 
we have also studied the correlation p^ z = p(log L* x , z). We note a positive coupling, pi z = 
0.53, for the full long GRBs sample, but - contrary to the ppT distribution discussed above - 
the piz values don't have an increasing trend, fluctuating between 0.36 and 0.61 for smaller 
u subsamples. For our upper envelope sample p^ z = 0.55, but it is accompanied with a large 
random scatter of (log L* x ,z) points. These features does not support a clear significant 
redshift evolution of the GRB afterglow luminosity distribution, but the issue should be 
studied in more detail (Dainotti et al., in preparation). Let us also note that our limiting 
upper envelope subsample includes GRBs with redshifts reaching the maximum value of 
'only' 2.75, while the most distant GRB with z = 8.26 disappears from the analyzed sample 
after decreasing u below 0.25. 

On the other hand to verify if the observed effect of higher LT correlation in the upper 
envelope is not a systematic effect due to higher photon statistics in brightest afterglows we 
compared distribution of the observed afterglow fluxes, F x , for the upper envelope sample 
with respect to the full analyzed GRBs sample, including also the IC GRBs (Fig. [5]). We 
find that the upper envelope observed flux values span more than 3 orders of magnitude 



- 11 - 



and are mixed in this range with other points, in the upper part of the flux distribution. 
The IC GRBs have the same behavior as the long GRBs, but they are on average less 
luminous. Thus, the L X (T*) and T* small errors result preferentially from the smooth light 
curve shapes, allowing precise fitting to the considered 'canonical' shape, not due to a higher 
observed photon statistics. We note in the considered distribution that the upper envelope 
points are preferentially in the its lower part, with the observed systematics resulting from 
the LT anticorrelation (for these GRBs fluxes are measured at larger fitted T*), but possibly 
influenced also by the mentioned weak L — z correlation in the sample. 



4. Summary 



In the presented analysis we discovered that the afterglow light curves which are smooth 
and well fitted by the considered canonical model, belong to most luminous GRBs forming the 
well correlated upper part of the (log T*, log L* x ) distribution. The GRB cases with appearing 
flares or non-uniformities of the light curves exhibit a trend to have lower luminosities for any 
given T*. We also noted the possible correlation of the X-ray spectral index /3 a and the time 
T*, which, together with the LT correlation, provide new constraints for the physical model 
of the GRB explosion mechanism. Let us also note that the revealed tight LT correlation, 
i f supported with lar ger statistics, could be a basis for a new independent cosmological test 
JCardone et alil2009h . 



An LT correlation for the independently analyzed (small) sample of IC GRBs is also 
revealed. It is different from the long GRBs, with a higher inclination of the fitted correlation 
lines and its luminosity normalization below the one for long GRBs. It provides a new 
argument for a separate physical reality of the postulated IC GRB sub-class. Any future 
attempt to study relations between various GRB properties should involve, in our opinion, a 
separation of the IC GRBs from the long ones, to limit analysis to physically homogeneous 
sub-samples (like our upper envelope one). A simple increasing the studied GRB sample 
with a mixed content may smooth out any existing relation. 

We do not intend to discuss here consequences of these findings for GRB physical models- 



Let us sim ply note that the LT r elation is predicted by the models of ICannizzo fc Gehrels 



(120091) and iGhisellini et alj (120091) . proposed for the physical GRB evolution in the time T a . 
The ICannizzo fc Gehrels I (120091 ) model predicts a steeper correlation slope (3/2) than the 
observed one (~ 1), which on the other hand is in a good agreement with the model of 



Yamazakil (120091 ) 
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Fig. 5. — L* x versus the measured energy fluxes Fx for the whole GRB sample. The 8 our 
upper envelope points are presented as red squares, while the IC GRBs are green triangles. 



